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We investigate the low-energy scattering and bound states of two two-component fermionic atoms 
in pure two-dimensional (2D) and quasi-2D confinements with Rashba spin-orbit coupling (SOC). 
We find that the SOC qualitatively changes the behavior of the 2D scattering amplitude in the 
low-energy limit. For quasi-2D systems we obtain the analytic expression for the effective-2D scat- 
tering amplitude and the algebraic equations for the two-atom bound state energy. We further 
derive the effective 2D contact potential for the interaction between ultracold atoms in the quasi-2D 
confinement. This effective potential can be used in the research on many-body physics of quasi-2D 
ultracold fermi gases with Rashba SOC. 
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I. INTRODUCTION 

The discussion on synthetic gauge field [l| and spin- 
orbit couphng (SOC) [3 in bosonic and fermionic sys- 
tems has recently drawn great attention [3| since its ex- 
perimental realization in cold atomic gases. By applying 
counter propagating Raman pulses with tunable proper- 
ties, the effective gauge field and SOC have been accom- 
plished in ultra-cold gases of both bosons @, Ij-Q and 
fermions jlG . Ill| . These experimental achievements add 
an additional piece to the already versatile toolbox of 
manipulation of cold atoms, and may help paving a new 
route toward the realization and investigation on novel 
quantum states. A considerable amount of theoretical 
interest has been stimulated in the understandin g of the 
SOC effect in both bosonic [H-i^ and fermionic [IS-lll 
systems. 

Among these works, one important direction is the 
study of SOC effect in low dimensionality where many 
interesting novel quantum states may be present. In par- 
ticular, a BEC with a half-quantum-angular-momentum 
vortex may exist in two-dimensional (2D) bosonic sys- 
tem with Rashba SOC [l7Hl9f . In two-component Fermi 
gases, a topological superfluid is proposed in 2D config- 
uration, and can supports zero-energy Majorana modes 
which are related to fault tolerant quantum computa- 
tion Hi, [13. 

In realistic experiments of ultra-cold atoms, the low- 
dimensional physics are usually studied in a quasi-low- 
dimensional configuration, where atoms are strongly con- 
fined in one or two spatial dimensions such that the trap- 
ping frequencies along these directions dominate all other 
relevant energy scales in an effective low-dimensional 
Hamiltonian. In the energy scale of interest, this effec- 
tive Hamiltonian should catch the same physics as the 
original Hamiltonian, which usually corresponds to two- 
body processes in the context of cold atoms as the par- 
ticle separation is much larger than the range of inter- 
atomic interaction. Thus, to write down the correct ef- 



fective low-dimensional Hamiltonian, one needs to inves- 
tigate the two-body physics in the quasi-low-dimcnsional 
confinement, and express the effective low-dimensional 
scattering amplitude or the two-body binding energy in 
terms of the "control parameters" , including the three- 
dimensional (3D) scattering length and the intensity of 
the confinement. In the absence of SOC, this analysis has 
been performed for quasi-one-dimensional (ID) [3^, 123] 
and quasi-2D j37l . |38| configurations, and leads to the 
unique feature of confinement- induced resonance [36[ and 
a different renormalization scheme around Feshbach res- 
onances [33, |33- In the presence of SOC, a discussion of 
two-body physics within quasi-low-dimensional confine- 
ment and a derivation of the effective low-dimensional 
Hamiltonian is still lack. 

In this paper, we investigate two-body scattering pro- 
cess and bound state energy of two-component fermionic 
atoms in 2D and quasi-2D geometries with a Rashba type 
SOC. For 2D systems, we find that when the total mo- 
mentum of two atoms is zero, the 2D inter-atomic scat- 
tering amplitude in the low-energy limit will be qual- 
itatively altered by the presence of SOC, rendering a 
polynomial rather than logarithmical decay to zero. For 
quasi-2D systems with a strong harmonic confinement 
along the axial z-dircction, we obtain analytic expres- 
sion for the effective 2D scattering amplitude. We also 
derive the equation for the binding energy of two-atom 
bound states (dimer) , as well as the effective mass of the 
dimer. We find that the presence of Rashba SOC tends 
to enhance the two-body binding energy. This observa- 
tion can be qualitatively understood by noticing that the 
density of states in the low energy limit is increased as 
the ground state becomes degenerate with SOC. Based 
on these results, we further derive an effective 2D inter- 
action between atoms in the quasi-2D gas with Rashba 
SOC, and map out two effective 2D Hamiltonians which 
are responsive for energy regimes around the two-body 
binding energy and close to the single-particle threshold. 
These effective models can be used to study many-body 



properties of quasi-2D gas with attractive or repulsive in- 
teractions, respectively. Our results provide the possibil- 
ity to control the effective 2D physics via 3D parameters 
including the scattering length and the intensity of z- 
confinement. The method developed in this manuscript 
can be directly generalized to other systems with all kinds 
of atoms and arbitrary types of SOC. 

The remainder of this manuscript is organized as fol- 
lows. In Sec. II, we discuss the two-body physics in 
2D. In Sec. Ill, we calculate the low-energy scattering 
amplitude and the binding energy of dimers in a quasi- 
2D geometry with a strong harmonic trapping potential 
along the axial z-direction. By matching the two-body 
physics, we construct in Sec. IV the effective 2D Hamil- 
tonian which can describe the low-energy behavior of the 
quasi-2D gas. The main results are discussed and sum- 
marized in Sec. V, while some details of our calculations 
are explained in the appendixes. 



II. 2D SCATTERING WITH RASHBA SOC 

In this section we discuss the scattering and bound 
states of two spin-1/2 fermionic atoms in a pure-2D ge- 
ometry [x-y plane) with Rasliba SOC. In such a system, 
the total momentum q of the two atoms is conserved. 
Then the spatial motion of the mass center of the two 
atoms is separated from the relative motion and the spin 
of the two atoms. Thus, the quantum state of the two- 
atom relative motion can be described by a spinor wave 
function 

IV (p)) = Vn(p)lt>ilt>2 + V^n(p)lt)iU>2 

+^it(p)U>ilt)2 + V'u(p)U)iU)2 (1) 

with p — {x, y) the 2D relative coordinate of the two 
atoms and | t)i(2) and | i)i(2) the spin eigen-states of the 
1st (2nd) atom. It is apparently that the wave function 
\^ (p)) can also be considered as a p-dependent spin stae 
of the two atoms. In this paper, we use the nature unit 
h = 7n = \ with m the single-atom mass. 

Due to the presence of SOC, the relative motion of the 
two atoms depends on their center-of-mass momentum 
q. For a given value of q ^ {QxTQy)i the Hamiltonian 
for atomic relative motion and spin degrees of freedom is 
given by 



^(2D)^^(2D)^^^^(^^^ 



(2) 



r(2D) 



H, 



(2D) 




Px 



where the 2D free Hamiltonian H^ ' is given by 

12 ^ 

3 = 1.2 

J = l,2 



'Y —+-y [-+(-1) 

2^ dp2^ 2 ^ 12 ^ ' 

I3=x,y 



rU) 



(3) 



with p = {px,Py) the relative momentum of the two 
atoms. The Rashba SOC is described by the second 



term of Eq. (|3]), where q/2 + (— l)-'p is the momen- 
tum of the j-th atom. The spin operator ax y for the 
j-th atom is defined as ax^ = | t)j(l I + I t)j(l I and 

ay^ = -i| t)j(l I + i| l)i(t I- Without loss of generality, 
we assume the SOC intensity ^ in Eq. ([3|) is a positive 
real number. 

In Eq. ^ V2D (p) is the atom-atom interaction poten- 
tial with p = \p\. We assume the effective range of the 
potential V2D (p) is p* , that is the potential becomes neg- 



ligible in the region p 



> 



We further assume the SOC 



is weak enough so that ^ ^ 4/p^. 

To investigate the scattering process in our systems, 
we first define single-atom spin state \aj,t)j for the j-th 
atom as 



? 



±(t,ay)+tya(f))\c 



• t),=a,^\a„t), (4) 



with t = {tx,ty) any 2D vector and aj = ±1. We further 
define the two-atom spin state \a {q, k)) as 

|a(q,fc)) = |ai,- + fe)i|a2,--fc)2, (5) 

with a = (011,02) and a = (a2,ai)- 

In the scattering process, the incident wave function 
shoud be the eigen-statc of the Hamiltonian H^ for 
the free motion of the two fermionic atoms. The straight- 
forward calculation shows that the such a incident wave 
function takes the form 



I^^"Up)) = 



23/27r 



\a{q,k))- 



-ikp 



23/27r 



\aiq-k)). (6) 



with incident momentum k ^ {kx,ky). In this paper we 
denote 



{a.,q,k) 



(7) 



as the set of all the three quantum numbers. Here we 
have considered the Pauli's principle for the fermonic 
atoms. The eigen-energy of H^ ' with respect to the 
cigen-state ji/'c (p)) is 



e 



"1 



k 



0L2 



k 



(8) 



with k = \k\. Notice that in the presence of SOC, the 
scattering threshold, or the minimum value of Sc with 
respect to a fixed q, is shifted from to £thrc (9) which is 
given by 



£thrc (q) = 



-q^i-e/^ {q<0 

-q^/2 (q > 



(9) 



with q ^ \q\. 

Now we consider the scattering state \ipc (p)) with 
respect to the incident state \^c (p))- We assume the 
scattering energy Ec is low enough with k <^ 1/p*- In 



such a low-energy case and within the region of p > p, , 

the wave function of the scattering state \ipc (p)) can 
be expressed as (see discussion in Appendix A and in 
Rcf. [H) 

|^i+)(p))«|V4"^(p)) + A(c)5(ec;p,0)|0,0), (10) 
where 

|o,o) = i/V2(|t)iU)2-U)i|t)2) (11) 

is the spin singlet state, and the 2D free Green's function 
9 im P, P') is given by 



9{mp,p') 



rj + iO+-H'o 



(2D 



:'5(P-P'), (12) 



and can be considered as a (p, p')-dependent operator for 
the two-atom spin. 

The coefficient A (c) in Eq. ([T0| can be derived with the 
following analysis. First, it can be proved (see Appendix 
B) that in the small- distance region p* ^ p -^ 1/fc the 
function |-0('^'+) (p)) behaves as 



|i/'^'''+^ (p)) oc (lnp-lnd)|0,0) 



(13) 



Here, the p-indcpendcnt factor Ind is determined by ^ 
and the detail of the potential V2D {p) (Appendix B). 
and almost independent on the scattering energy Ec in 
the low-energy case [38|. Second, the calculations in Ap- 
pendix C shows that, in the small-distance region we have 



(0,0|.g(r7;p,O)|0,0) « ^ 



Inp+C+ln ( -z-^-- 



^A(ec,q') 



(14) 



where C — 0.5772... is the Euler gamma number and the 
function A(r/, q) is defined as 



A(^, 



(27r) , J 



?; + iO+-ec" 77 + iO+-|fc 



"12 



(15) 



with c" — {a.";q;k"). For q = 0, a straightforward 
calculation shows that the function X{ri, 0) can be further 
simplified as 



Hv,o) 



2(27r) 



2^/-iJl-e 



TT + 2 arctan 



C+2v 



in < -^) 



_i^5 _^^arctanh(2V'?+eV4/c) 



2aA+?V4 



Vv+Wi 



(-^<77<0) 



—pd arctanh — ^J : (0 < 77 < 1 - 4r) 



(16) 



Substituting Eq. (fT4|) into (TTOl) . we get the expression 
for I -0c (p)) in the small-distance region 

IV'^) (p)) « iv4"^ (p)> + 
41r{i"^+^+i"(-¥)+^["-9]}|0'0)^i7) 

Comparing Eqs. (|17p and (J13p . we obtain the expression 
for the parameter A (c) : 



A(c) 



(2^)(0,0|4''^(O)) 



iTT/2 - C - In (d^/2) - (27r)A (ec, g) ' 



(18) 



According to Eq. (|10p . A (c) completely determines the 
behavior of the scattering-state wave function |f/;('^'+) (p)) 
in the region of p > p* . 

According to scattering theory [43], we can define the 
2D scattering amplitude Z'-^'^-' between the incident state 
\ipc (p)) and an energy-conserved output state \ipl., ' (p)) 
with q' ~ q and Ec = Ec' as 

/(2D) (,' ^ ,) = _2^2 /dp(^,W (p) |1^2d|0(+) (p)). 



(19) 
A straightforward calculation (Appendix D) shows that 
we can express f^"^^^ in terms of the coefficient A (c) 



/(2D) (c'^c) = -2^2^ V'^?^(O)|0,0) A (c) 



(20) 



Now we discuss the two-body bound states with 
Rashba SOC in two-dimensions, as is also investigated 
in Ref. |34| . As shown in Appendix A, when the energy 
£(, of the bound state is close enough to the scattering 
threshold etiirc('z), or the condition ethrc((z) — £& ^ 1/p* 
is satisfied, the wave function j-i/jf, (p)) of the two-atom 
bound state can be approximated as 



|7Ab(p))«B.g(e6;p,O)|0,G) 



(21) 



in the region p ^ p*. Here, B is the normalization coef- 
ficient and g is defined in Eq. (|12p . The energy ei, of the 
bound state is determined by the condition (appendix B) 



(0,0|V'6(p)) oclnp-lnd 



(22) 



in the region p* < p <^ ^/Vl^b - ethrc (g)|, or by the 
equation 



-lnd=C + ln - 



V^fc 



(27r)A(£b,q). 



(23) 



Here, we use the fact that in the small-distance region the 
function (0, 0|(7 (£&; p, 0) |0, 0) also takes the form as in 
Eq. ([13]), with £c replaced by the new variable Sb- There- 
fore, the bound-state energy Sb is a function of both the 
characteristic length d and the center-of-mass momentum 

In the discussion above, we obtain the analytical ex- 
pressions for the scattering amplitude f^'^^' and the equa- 
tion for the bound-state energy in a pure 2D geometry 




e(c)(units of d ) 



FIG. 1: (Color online) The variation of function |_Fp with 2D 
scattering energy e{c). Here, F is defined in Eq. (|26[l and d is 
the parameter in Eq. p3p . We show results in cases with zero 
center-of-mass momentum q = and SOC intensity ^ = (red 
solid line), 1/d (blue dashed-dotted line), 2/d (black dashed 
line) and 5/d (green open circle). 



We emphasize that, this change is due to the dispersive 
relation (|5|) of Ec- When q = 0, Ec is independent of the 
direction of fc, and takes minimum value — ^^/4 when 
k = ^/2. Namely, we have Ec = — C^/4 for all momenta 
fe in the cirecle of {kx,ky) plane. Nevertheless, when 
q 7^ 0, or when SOC is not of Rashba type, this character 
disappears. In these cases Ec takes minimum value only 
when fc takes one or two certain values, as in the systems 
without SOC. Thus, the 2D scattering amplitude also 
logarithmically decays to zero when the scattering energy 
approaches to the threshold. 

Second, when there is no SOC, the scattering ampli- 
tude is independent on the center-of-mass momentum q. 
This feature is qualitatively altered by the presence of 
SOC, as can be clearly seen from Eq. ^ where the SOC 
terms couple the relative motion of the two atoms with 
the total momentum q. As a consequence, the scattering 
amplitude /'■^^■' becomes a function of q. The similar 
q-dependence can also be observed in the bound-state 
energy Et of the low-energy bound states [3J] . 



with Rashba SOC. Comparing our results with the 2D 
scattering theory without SOC (see Appendix E and Ref. 
^^), we observe the following two qualitative differences: 
First, when the total momentum q of the two atoms 
is zero, the SOC changes the dependence of the scatter- 
ing amplitudes on the scattering energy Ec- As shown in 
Appendix E, when there is no SOC, the scattering am- 
plitude /g decays to zero logarithmically in the limit 



0, 



liin /( 



(2D) 



cx 



ln£. 



(24) 



When q = 0, a Rashba SOC will change the scatter- 
ing amplitude through the function A (ec, 0) in A (c) and 
the factor (V'^?^ (0) |0, 0)(0, 0|?a1^°^ (0)). In particular, the 
A-function removes the logarithmic behavior of the scat- 
tering amplitude in the region around Ec — 0, leading 
to 



lim f 



(2D) 



OC 



4 ■ 



(25) 



III. SCATTERING IN A QUASI-2D 
CONFINEMENT WITH RASHBA SOC 

In the previous section, we obtain the two-atom scat- 
tering amplitude and bound-state energy in a pure 2D 
system with Rashba SOC. Our results show that the SOC 
qualitatively changes the 2D scattering amplitude. In a 
realistic experiment of cold atoms, the 2D condition is 
usually realized by applying a tight confinement along 
one (say z) direction, such that the degrees of freedom of 
the single-body Hamiltonian along that specific direction 
can be approximately integrated out. In this section, we 
study the scattering and bound states of two spin-1/2 
fermionic atoms in a quasi-2D configuration. We con- 
clude that the qualitative effects introduced by Rashba 
SOC to the two-body physics in pure 2D geometry also 
exist in the quasi-2D cases. Besides, we show how the fac- 
tor In d in the expressions of 2D scattering amplitude and 
bound state energy can be effectively tuned by the atomic 
3D scattering length and the intensity of z-confinement. 



Namely, in the presence of Rashba SOC, the scattering 
amplitude f^^^' polynomially decays to zero in the limit 

of Ec ^ -e/4- 

To illustrate the SOC effect to the 2D scattering am- 
plitude, we plot in Fig. 1 the mode square of the quantity 



F : 



/(2D) (c' ^ c) 



,,(0) 



{i;]::>{o)\o,o){o,om"m 



(0) 



(26) 



as a function of the scattering energy Ec with respect 
to different SOC intensities. It can be seen clearly that 
the function |_F| linearly decays to zero in the limit of 
Ec — )■ — C^/4 in the presence of SOC, and logarithmically 
decays to zero when there is no SOC. 



A. System and Hamiltonian 

We consider a quasi-2D configuration with a harmonic 
trap with frequency ui applied along the z-direction, while 
the atomic motion in the x-y plane is free. In the quasi- 
2D case, the total momentum q of the two atoms in the 
x-y plane is also conserved and serves as a parameter for 
the atomic relative motion, as in the pure 2D case. For a 
given value of q, the Hamiltonian of the atomic relative 
motion and spin states takes the form 



H = H, 



(2D) 



H, + 1/3D {r) . 



(27) 



Here, Hq ' is defined in Eq. ((3)), 



H, 






"2 



(28) 



is the Haniiltonian for the two-atom relative motion in 
the z-direction, with z the relative coordinate of the two 
atoms in this direction. Note that we have shifted the 
zero-energy point for our convenience, since the harmonic 
confinement allows us to separate the relative motion in 
the z-direction from the center-of-mass degree of free- 
dom. 

In Eq. (j27p . the atom-atom interaction potential 
V3D {r) is a function of the relative position r = (p, z) = 
{x,y, z) of the two atoms in three dimensions. Here, we 
have r = |r| and denote the effective range of ^30 (r) as 
r^,. In the region r >^ r^,, the potential V3D ('') becomes 
negligible. For simplicity, we further assume V3D [r) is 
independent on the atomic spin, and consider only the 
weak SOC case with 4/^ » r*. 



B. Two-atom scattering state 

In the quasi-2D configuration, we consider only the 
scattering between two fermionic atoms in the ground 
state of Hz ■ Then the incident wave function is given by 

|vI/(o)(r))=^oWIV4"^(p)> (29) 

with \ipc (p)) defined in Eq. ([5]). Here, (pn, (z) is the 
Uz-th eigen-wave function of Hz ■ Without loss of gener- 
ality, we set the phase of (/3„^ (z) so that fn^ (0) is real. 
We further assume the energy gap between the incident 
state and the threshold ethro is smaller than the trapping 
frequency along the z-direction, i.e. \ec — ethro('7)| < w. 
In this case, the energy-conserved output states of the 
scattering are also in the ground transverse channel with 

Next, we calculate the scattering state |\l/c (r)) with 
respect to the incident state I'i'c (r)). The scattering 
state can be obtained with the same method as in Sec. 
II. As shown in appendix A and Ref. [Sa], in the region 
of r > r*, the scattering state wave function |^c (r)) 



|vl/(+)(r.))«|vI/(0)(r)) 



Acs (c) 
'/'o(O) 



G(£e;r, 0)|0,0), (30) 



with the quasi-2D free Green's function G (??; p, p') given 
by 

1 



G{mr,r') 



T] + iO+ - 



H 



(2D) 



Hz 



6ir-r'). (31) 



Here A^s (c) is a r-independent coefficient, and can be 
derived with the following two facts. First, it can be 
proved (see Appendix B) that in the small-distance re- 
gion r^, < r <$i 1/fc the function j'l'c (r)) behaves as 



(0,0|*(+)(r)) 



1 1 

oc I 

r a 



(32) 



with a the s-wave scattering length determined by ^ and 
the detail of V3D {r) (Appendix B). Second, Appendix C 
also shows that, in the small-distance region we have 



(0,0|G(77;r,O)|0,0) 



1 

Anr 



00 

X! \'PnMfH£c-nzU};q), 



(33) 



n-=0 



where l^ = -^/l/w, w (ry) is defined as 
w{r]) = 



N N 
lim I 2\ — In ^r 

N^oo \ V TT e^ 



N 



^Rizill! ,„(,-,-.„. 



j=o 



(2j)! 



and A(?7; q) is defined in Eq. (fT5|) . Substituting Eqs. ([30|) 
and ([551) into ([5^ . we can obtain the parameter A^ff (c). 
We find that A^s (c) can be formally expressed as in 
Eq. (fT8)) . with the 2D effective range d replaced by a 
quasi-2D effective range doff 



Acs (c) 



(2^)(0,0|4"^O)) 



/2 -C -In {des (ec, q) V^c/2} - (27r) A (e^, q) ' 



ITT 



(35) 



Here, the effective characteristic length d^s is a function 
of the scattering energy Sc, the center-of-mass momen- 
tum q, and the 3D scattering length a, and can be de- 
termined by 



Indoff (ec,<7) = - 



2™(ft; i,f_H^|_c 



2 



-^ + {2nf lo V |^„, (O)l'A(e, - n.c.; q). (36) 

Note that the summation of Uz on the right-hand-side of 
Eq. p6p runs over all natural numbers. One can easily 
show that Indcs {r]-,q) always takes a real value when 
1] < etiirc(9) + w. We would like to emphasize that the 
term In d^s governs the effective 2D physics in this quasi- 
2D system, and includes all effects given by the control 
parameters a and uj. 



C. Effective 2D scattering amplitude 

The straightforward calculation shows that, in the re- 
gion p ^ Iq wc have 

|vl'(+) (r)) « ^0 (z) (|^(°) (p)) + Aeff (c) g (£,; p, 0) |0, 0)) , 

(37) 
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FIG. 2: (Color online) The variation of l-Foffl^ as functions 
of the quasi-2D scattering energy e(c). Here, F^a is defined 
in Eq. (|40[) and /o = l/y/ZJ. We sliow the results with zero 
center-of-mass momentuni g = 0, SOC intensities ^ = (a) 0, 
(b) 0.3/^0, (c) 1/^0, (d) 3//o and 3D scattering lengths a = oo 
(red solid line), lo (black dashed line) and —Iq (blue dashed- 
dotted line). 
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FIG. 3: (Color online) The variation of |Feff|^ with 3D 
scattering length. The function F^s is defined in Eq. (|40p 
and lo = Xj'JbJ. We show the results with zero center- 
of-mass momentum q = 0, SOC intensities ^ = (a) 0, (b) 
0.3/Zo, (d) 1/^0, (d) 3/^0 and the quasi-2D scattering energy 
£(c) = -CV4 + 0.006tJ (red sohd line), -^ j^ + 0.02tj (black 
dashed line) and — ^^/4 -I- 0.2a; (blue dashed-dotted line). 



where the 2D Green's function g is defined in Eq. ([T2|) . 
Comparing Eq. ((37|) and Eq. (fTO|) . we find that in the 
long-range region with p '^ Iq, the quasi-2D scattering 

state wave function |$c (j")) is the product of fo{z) 
and a 2D scattering-state wave function. Therefore, if 
we focus on the long-range region, the quasi-2D scatter- 
ing process is equivalent to a 2D scattering process with 
an effective characteristic length doff (^c, q), which is con- 
trolled by the 3D scattering length a and the trapping 
frequency a; in the z-direction. 

The quasi-2D scattering amplitude between the in- 
cident state l^c (r)) and an energy-conserved output 
state 1^^/ (t)) with q = q' and Ecj — £c is defined as 

/(2D) (c' ^ c) =-2^2|^r'(*(?) (r') | V^r') \¥+^ (r')). 

(38) 
As shown in appendix D, /(Q^^) (c' <— c) can be ex- 
pressed as 



/(Q2D) (^, ^ ^) ^ -2^2^7^(0) (Q)|o^ Q^^^^ (^) 



(39) 



Comparing Eq. 1^ with ^, we find that /(Q2D) jg 
nothing but the 2D scattering amplitude with respect to 
the 2D scattering state given by the right-hand side of 
Eq. dSll). 

To understand the behavior of the effective quasi-2D 
scattering amplitude /(Q^d)^ ^^ show in Figs. 2 and 3 
the variation of the mode square of the function 



F, 



cff 



j(Q2D) (^, ^ ^) 



.,(0) 



(V'i;'^(o)|o,o)(o,oiV'r(o)) 



(0), 



(40) 



with scattering energy Sc, SOC intensity ^, 3D scattering 
length a and characteristic length Iq of the z-confinement. 
In Fig. 2, it is shown clearly that the quantity |i^cff | lin- 
early decays to zero in the low-energy limit £c — >■ ethro('z) 
with SOC, and logarithmically decays to zero when there 
is no SOC. This observation is consistent with the out- 
come in the pure 2D systems as discussed in the previous 
section. In Fig. 3 we investigate the behavior of |i^off| as 
functions of the scattering length for different values of 
£c — £thro(9)- Note that for a given value of £c — £thic((z), 
the resonance behavior of |i^off | is still maintained, while 
the resonance point is shifted by the SOC and the am- 
plitude of |Foff| is suppressed by the SOC. 



D. T"wo-atom bound state 

Next, we consider the two-body bound state in quasi- 
2D configuration. Similar as in Sec. II, it can be shown 
that in the region r > r,, the wave function |\l/b (r)) can 
be expressed as 



\^b{r)}=B'G{Eb;r,O)\0,0) 



(41) 



with B' the normalization factor. In the long-range limit 
r ^ lo, the wave function of the quasi-2D bound state 
is proportional to that of the pure-2D bound state in 
Eq. (PT|) . The energy Et of the quasi-2D bound state 
is determined by the boundary condition in the small- 
distance region (appendix B) 



(0,0|*b(r))cx 



1 
a 



(42) 
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FIG. 4: (Color online) The binding energy -Ebinding = 
ethrc(g) — El, of the quasi-2D two-atom bound state as a 
function of the 3D scattering length. In this plot, we set 
the center-of-mass momentum q — and the SOC intensi- 
ties £, — (black dots), 0.5//o (blue dashed line), l/lo (green 
dashed-dotted line) and 3/lo (red solid line). 



Then it is easy to prove that the boundary condition 
Eq. (P^ is equivalent to the equation 



\nd^s{Eb,q) = C + ln - 



27:X{Eb,q). (43) 



In Fig. 4, we show the behavior of binding energy 

-^binding = £thro('Z) — Eb (44) 

with varying 3D scattering length a and SOC intensity 
for cases of g = 0. Notice that the variation of -^binding 
in terms of 1/a has the same quahtative behavior with or 
without SOC, and the value of i^binding is increased with 
the SOC intensity ^. 

Next, we discuss the dispersion relation of the two- 
atom bound state. To this end, we express Ei, as 
Eb = X)^o Ebnq"hn (q/q) and then substitute this ex- 
pression into Eq. ((33]). Expanding both sides of Eq. (pS]). 
we can obtain all the coefficients Ebn and the functions 
hn- When q is small, we have 



Eb ~ Ebo + Eb2q 
with EbQ and Eb2 determined by the equations 



(45) 



IoV2tt 
2a 



(2^) 



1 



EbO^ 

"2^ 



Vn, (0)1 Ao(£'bo - n^uj) 



(46) 



and 



E, 



62 



J2Z=f)\fn. {0)\^>^2{EbO - n^uj) 



'^'{^hEZ^olVnMrKiEbo-n.^) 



4(2ir)3/2c^;o 



(47) 




3.5 



2.5 



(b) 
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FIG. 5: (Color online) (a) The variation of the effective mass 
m,cff of the two-atom bound state in the quasi-2D system as 
a function of the 3D scattering length. The SOC intensities 
used in this plot are ^ = l/lo (black dashed line), 2/lo (blue 
dashed-dotted line) and 5/lo (red solid line), (b) The variation 
of TTicfi as a function of the SOC intensity ^. The 3D scattering 
lengths used in this plot are a — —lo (black dashed line), oo 
(blue dashed-dotted line) and lo (red solid line). 



Here, the functions Ao (??) and X2{ri) are defined as 



and 



^oiv) 



A2(^) 



M TT + 2 arctan 



e+2v 



8ny/-4T]-£_^ 



(48) 



^ ^^-477-^2 _^_ 2(2,; + ^2) arctan 



V-4')-?" 



32^77 (-477-^2)^^' 



(49) 



with w' (77) = dw (77) /d77 and Aq(77) = dAo (7;) 7^77. 

The total energy of the two-body bound state can then 
be obtained by adding Eq. (|45|) and the kinetic energy of 
the centcr-of-mass motion, leading to Ebt ~ EbQ+Eb2q'^ + 
q^ /A. This quantity can also be expressed in terms of the 
effective mass 



1 



Weff 



2E, 



62 



1/2' 



and takes the form 



Eb2 ~ Ebo 



2mcff 



(50) 



(51) 



In Fig. 5, we plot the effective mass with the 3D scatter- 
ing length and the SOC intensity. 



IV. EFFECTIVE 2D HAMILTONIAN 

With the above knowledge of two-body physics, we 
can construct the effective 2D Hamiltonian for a two- 
component Fermi gas with Rashba SOC in a quasi-2D 
confinement. The effective Hamiltonian is required to 



give the same scattering amplitude or two-body bound 
state as the original Hamiltonian. It is pointed out that, 
for each effective 2D Hamiltonian, this criteria can only 
be satisfied for a small energy range. Thus, in principle 
one needs to derive different effective Hamiltonian for 
different energy regions. 

We first consider gases with atoms in the low-energy 
scattering states i.e., the system given by the directly 
cooling of the atoms with a fixed scattering length. In 
such a system, the energy of relative motion of two atoms 
is slightly above — ^^/4, and the probability for atoms in 
the excited states along the z-direction is negligible. The 
atom-atom interaction can be described by a 2D single- 
channel contact potential 

^o = ^Yl '5(fc + fc')aLt4',i«fc'+fc"4«fc'-fc",f (52) 

k,k' ,k 

Here, a|, ^ and Op o- are the 2D creation and annihilation 
operators for an atom with momentum p and spin a, S 
is the area of the system, and the summation ^^ ^; j^.// is 
done for |fc - k'\/2 < k^ and |(fc - k')/2 + k"\ < L, with 
kc a cut-off momentum. For systems without SOC, the 
coupling intensity g{q) is given by the renormalization re- 
lation ^ l/(27r.g) = -Xc<j. dk{Eb + Py^, with Eb the 
bound-state energy. For our current systems with Rashba 
SOC, the renormalization relation should be modified as 



Infer — 



2^g{q) 



-C-ln 



4ff -i-'9 



(53) 



with the effective characteristic length d^s defined in 
Eq. p6)) . According to the definition, doff also depends 
on 3D scattering length a and the characteristic length 
Iq of the z confinement. A straightforward calculation 
(Appendix F) shows that the 2D scattering amplitude 
given by the interaction potential Vo is the same as the 
quasi-2D scattering amplitude /(Q^d) ^^ -^^ ^^ ([M|) . 

Now we consider the systems with atoms in the bound 
states. In that case, if the binding energy of the bound 
state is large, the atomic population in the transverse 
excited states <fn, [z) with n^ > becomes significant |37l . 
l4l| . To take into account these transverse excitations, we 
use a phcnomcnological two-channel model to describe 
the atom-atom interaction: 

q ^ ' 

+ ^= ^ 'u{k + k!)a\ ^a\, ^bk+k' + h.c. (54) 
^^ kM 

Here, the open-channel interaction Vo is defined as in 
Eq. (j52[) . 6]- and bq are the creation and annihilation 
operators for molecules in the closed channel of our phc- 
nomcnological model with [6q,6L] = Sq,q', and the sum- 
mation X)fc fc' is done for |fc — k'\/2 < kc- The renor- 
malization relation for cases without SOC is given in 



Ref. [4l|. In the presence of Rashba SOC, The parame- 
ters v(q) and u(q) are now given by 



1 



K{Eb,q) 



.9(9) 



uiqy 



Eb - v{qf 



(55) 



and 



Ei,-v(q) 



5(g) 



<qy 
Eb~v{q) 



j:Z=i\^r.AO)\ xiEb-n,L,,q) 



(56) 



where the functions k {rj, q) and x iVi q) are defined as 



K (77, q) = -27r In fee + TT In |?7| -f (27r)^ A {?], q) (57) 



and 



X iv, q) 



E 



,^\m^i^^ (58) 

k<k, (77 - Ec) 



with c = (a, q, k). In Appendix F, we show that the 2D 

bound state 1^1'^^° (t*)) given by Vb has the same binding 
energy as the quasi-2D bound state |^b(r)) obtained in 
the previous section. Besides, it is also proven that the 
open-channel probability of \^^ (r)) recovers the popu- 
lation of the transversal ground-state |*;,(r)). Therefore, 
the effective 2D potential VJ, can be used to study low- 
energy physics where nearly all particles are in the bound 
states with zero ccnter-of-mass momentum. 



V. CONCLUSIONS 

In this paper we investigate the two-body physics of 
two-component fermonic atoms in 2D and quasi-2D con- 
figurations with a Rashba SOC. For the 2D case, we 
find that when the total momentum is zero, the loga- 
rithmic behavior of the 2D scattering amplitude in the 
low-energy limit is replaced by a polynomial behavior. 
For the quasi-2D system, we obtain an analytic expres- 
sion of the effective 2D scattering amplitude as a function 
of the 3D scattering length a and the trapping frequency 
Lu along the strongly confined z-direction, and observe the 
same polynomial behavior. We also discuss the two-atom 
bound state, and get the algebraic equation for the bind- 
ing energy. We find that the two-body binding energy 
is enhanced by the presence of SOC, as a consequence 
of the increase of single particle density of states in the 
low energy limit. We also analyze the dispersion relation 
and extract the effective mass of the dimers. These in- 
formation allows us to tunc the effective 2D physics in 
quasi-2D configuration with parameters a and w. 

With the knowledge of two-body physics, we further 
construct two effective 2D Hamiltonians, which can indi- 
vidually mimic the original quasi-2D Hamiltonian within 
the energy regimes around the two-body binding energy 
or close to the single-particle threshold. These effective 



models can be used to analyze many-body physics of the 
system, hence help paving a way towards the simulation 
of 2D Fermi system with quasi-2D atomic gases. Our 
method developed in this paper can be directly general- 
ized to systems of bosonic or distinguishable atoms with 
arbitrary type of SOC. 



Acknowledgments 

We thank Xiaoling Cui and Hui Zhai for helpful dis- 
cussion. This work is supported by National Natural 
Science Foundation of China (11074305, 10904172), the 
NKBRSF of China (Grant No. 2012CB922104), the Fun- 
damental Research Funds for the Central Universities, 
and the Research Funds of Renmin University of China 
(10XNL016). PZ and WZ would also like to thank the 
NCET Program for support. 



Appendix A: wave functions of 2D scattering states 
and bound states 



In this appendix we derive the Eqs. (JTUIIHU)) for 2D and 
quasi-2D low-energy scattering states, as well as Eqs. (|2n 
l¥T|) for 2D and quasi-2D bound states. 

We first prove Eq. pU)) for the 2D scattering wave func- 
tion \il)c [p)) ■ The Lippmann-Schwinger equation gives 

\^'t^ (P)) = l^f {p))+ldp'g{e,;p,p')VMp') \4'''> (p'))- 

(Al) 
Here, \ipc (p)) is the incident state and the free Green's 
hmction g{ec; p, p') is defined in Eq. (fT2|) . Since the po- 
tential V2D (p) is negligible in the region p ^ p*, the inte- 
gration in Eq. (jAip is only effective in the region p' ^ p* . 
On the other hand, in the low-energy cases k ^ l/p*, 
when p ^- 00 and p' < p*, the function g{ec]p,p') be- 
comes very steady with respect to p' and we have 



g{£c;p,p') «g(ec;p,o). 



(A2) 



Therefore, in the limit p -^ 00, the solution of Eq. (|A1[) 
takes the form 



l^i+'(p)) = IV^^"^(p))+5(ec;p,o)|x), 



(A3) 



/,(+) 



where the spin state |x) is related to \ipc (p)) via the 
equation 



Ix) 



dp'vMp')\rc''Hp')) 



(+) 



(A4) 



On the other hand, since \ipc (p)) is an eigen-statc of 



r{2D 



-^0 (9) + ^2D (p) and the potential V2D (p) is negligible 
in the region p ^ p*, in such a region the wave function 

\ipc (p)) satisfies the equation 



Therefore, the behavior of the wave function | V'c (p)) 
in the region p ^ p* is determined by Eq. (|A5[) and the 
boundary condition (|A3|) in the limit p — >■ 00. It is easy 
to prove that, the function |-0c ) + gi^c', P, 0)|x) satisfies 
both of the two conditions. Therefore, IV'c (p)) satisfies 
Eq. (jASP not only in the limit p — >■ 00, but also in the 
entire region of p > p* . 

Furthermore, due to the facts Pi2^2d(p) A2 = V2d(p) 

and Pi2\ipi'^^ (p)) = -IV'c^^ (p)) with A2 the permuta- 
tion operator of the two atoms, the integration (jA4p can 
be re-written as 



dp'v2B{p')\4''Hp')) 

= |0,0) fdp'V2u{p'){0,0H'tHp'))- (A6) 



(A7) 



Thus, Ix) can be expressed as 



\x)=A{c)\Q,0) 



withyl(c) a c- number. Substituting Eq. (jA7p into (|A3|) . 
we finally obtain Eq. ([10]). 

It is easy to find that Eqs. ^ ^ and (gl]) can be 
proved within the same approach. Especially, we have 
the result 

^^^|0,0) = fdr'VMr') l*(+^ (r-')>- (A8) 

V'o (0) J 

which is similar as the ones in Eqs. (|A4p and (|A7I) . 



Appendix B: small-distance behavior of 2D and 
quasi-2D wave function 

In this appendix we prove Eqs. dH El [H |12|) for the 
behaviors of the 2D and quasi-2D wave functions in the 
small-distance region p* << p << 1/fc or r* << r << 
1/k. Here we only show the proof of Eqs. ([T^I5^ . Eqs. 
(Pl 22) can be derived with the same method. 

We first prove Eq. p^ for the behavior of \ipc (p))- 
To this end, we define a rotated wave function 

|^(+)(p)) = f/(p)|V^W(p)) (Bl) 

with the unitary transformation U (p) defined as 



U (p) == cxp 






(4^) - -^^ 



exp 



i£,y 



r(l) 



Then the eigen-equation 

ij(2°)|V'W(p))=ec|V'i+)(p)) 
satisfied by |V'c (p)) gives 



V 

(B2) 



(B3) 



nn^i-^Hp)) 



sM^Hp)). 



(A5) 



ij(2D)|^(+)(p))^,^|^(-h)(p))^ 



(B4) 
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where 



H^^^^ = U{p)H'^^''^U^{p) 



^ 9/32 

P=x,y ' 



i^-^W (P) (^ ) + ^' (9' P) + ^2D (P) 



_ fV(2D) 



i?, 







V^2D (p) 



(B5) 



Here, we have V2D = UV2dW and the operators W^ (p) 
and W^' (<7, p) are defined as 



and 



(!)_,(.) 



M/'(q,p) 



16 ^^ 

a.—x,y 
a— j;,y 



+ C[/t(aW-a(2))[/ (B6) 
t/^(p)(^i^^--i^^)f/(p) 



16 



^|«'(.). 



(B7) 



Due to the weak SOC condition ^ << 1/p*, in the region 
p <C 4/^ we have we have W {p) ~ and VF' (q, p) ~ 
VK'(q, 0). Therefore, in this region the rotated wave 
function \ipc (p)) is determined by the Hamihonian 



E 

P=x,y 



a/32 



W (g, 0) + V2D (p) 



(B8) 



without SOC. Then the behavior of ji/Jc (p)) in the 
small-distance region is the same as the one for the 2D 
scattering wave function between two fcrmions without 
SOC, and can be described as 



|Vi^+^(p)) ex (lnp-lnrf)|0,0). 



(B9) 



with the characteristic length d determined by the 
detail of 1^d(p)- Using the relation j^/^c (p)) = 
U^ (p) |i/Jc (p)), we immediately get the result in Eq. 

(fT3| for the small-distance behavior of \ipc (p))- 

Eq. ([5^ can be proved with the similar approach. We 
can define the quasi-2D rotated wave function 

|*i+)(r)) = C/(r)|vI/(+)(r)). 
Then |^i+^ (r)) satisfies 



(BIO) 



H, 



(2D) 



+ //. + V^3D (r)J |*i+) (r)) = £,|^(+) (r)), 

(Bll) 

with #^^°' defined in Eq. (|B5| and ^30 = UV^^U'' . 
Similar as above, in the region r ^ 4/^ we have W (p) ~ 

and ly' (q, p) sa M^' (q, 0), and 1*^'*"' (r)) is determined 
by the Hamiltonian 

^ E ^ + H. + W'{q,0)^%T,{r) (B12) 



I3=x,y 



9/32 



without SOC. Then the behavior of \^i'^^ (r)) in the 
small-distance region is the same as the one for the 3D 
scattering wave function between two fermions without 
SOC, and can be described as 



|^W(r))(x(^i-i)|0,0), 



(B13) 



with the scattering length a determined by the de- 
tail of V3d(^)- Using the relation {"^Sli^^r)) = 
U^ (r) l^i^-* (r)), we obtain the behavior of l^i^-* (r)) in 
the small-distance region: 

|*WM)cxfi4)|0,0> + f E7(-?'--?0lO'«)■ 
(B14) 



l3=x,y 

Then wc have the result in Eq. (^^ 



(0,0|vI/^+)(^))^ 



1 1 



(B15) 



In the end of this appendix, we point out that, the results 
in Eqs. (jB9[) and (|B13[) of this appendix can be proved 
more explicitly with the approach given in Ref. [43| . 



Appendix C: small-distance behavior of 2D and 
quasi-2D Green's function 

In this appendix we prove Eqs. ((Til [55|) for 
the behavior of the functions (0,0|(7(ec; P, O)|0,0) and 
(0, 0|G(ec;T",O)|0, 0) in the small-distance region. We 
begin from Eq. ([H]). In the small-distance region 
P* ^ P « 1/^; the function (0, 0|g(£c; P, 0)|0, 0) is gov- 
erned by the leading terms in the limit p — >■ 0. Using the 
facts 

5{p-p') = j dk^--^^ f E I" (9' fc)> (" (9' fc) I ] 

(CI) 
and 

ij(^°) {e^'^P\cx [q, k))) - e, {e^'^P\cx (q, k))) , (C2) 

it is easy to show that 

(0,0|.g(£c;p,O)|0,0) 



E/*' 



r,ik' -p 



{2nY 



dk' 



Jk'p 



{0,0\a' {q,k'))\ 

Ec + iO+ - Ec' 
1 



(2^)'ec + ^0+-|fc'|2 

ik' p 



5^/dfe'^|(0,0|a(q,fc')>| 



iO^ 



tO^ 



(C3) 
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with c' = (a'; q,k'). Using the fact 



dk'- 



(C4) 



(27r)'£c + ?:0+-|fc'|2 277 

with Xo the modified Bcssel function, wc get the result 
(0,0|g(£e;p,O)|0,0) 
;A'o {-i^/e'cp) 



27r 



cik' p 



X /dfc'^|(0,0|a(q,fc')>| 



ec + iQ+-ec' £c + iO+-|fcT 



In the sniaU-distance region we have 



2 hip+C+ln , 2 



(C5) 



and 



Jk'p 



J2 dk'-—,\{OA^{q,k'))\ 



1 



1 



A(ec,9). 



(C7) 



with A- function defined in Eq. (|15p . Then we have proved 
Eq. dUl). 

Eq. psp can be proved with the similar approach. In 
the small-distance region r^, < r « 1/fc, the function 
(0, 0|G(ec; p, 0)|0, 0) is governed by the leading terms in 
the limit r ^ 0. Using the facts 



6{r~ r') 



dk 7^ — I y ^ \a (g, k)){a {q, k) 



{2^Y 



^ (y5„^ {z) ^Pn, {z) \ 
\ n J 



(C8) 



and 



= (nLc- + e,)(^„^(^)e''=-''|a(q,fc)), (C9) 

it is easy to prove that 



lim(0,0|G(ec;r,O)|0,0) 

r— fO 

1 

lim 



r-s-O 



Ec + iO^ 



2-^OL^x.y 






Sir) 



Yl IfnMfH^c-n.uJiq) 



(CIO) 



As shown in Ref. [38[, we have 
lim 



1 1 



iO^ 



2-^Qi- 



-x,y da'^ 



H, 



S{r) 



w{ej2) 

47r r 2 {2nf/^ Iq 



(Cll) 



with the function w (rj) defined in Eq. p4p . Substituting 
Eqs ((CTT|) into Eq. ([UTO)) . we can get Eq. (Pj). 



Appendix D: Scattering amplitudes in 2D and 
quasi-2D geometries with SOC 

In this appendix we proof Eqs. (PO)) and ([55)1 for the 
2D and quasi-2D scattering amplitude. The 2D scatter- 
ing amplitude Z'^^-* (c' <— c) between the incident state 
{ipc (p)) and out-put state |^^, ' (p)) is defined as 



(C6) /(2D) (,. ^ ,) = -2n^ldp{i,[?^ (p) I y2D(p) 1^-^+) (P)), 

(Dl) 



In the region p' < p, we have 



(D2) 



Substituting Eq. (|D2j) into (|D1[) . we immediately obtain 
the resuh in Eq. (pi)) : 

/(2D) (c' ^ c) = -27r2(v4?) (0) |0, 0) A (c) . (D3) 

Eq. ((39)) for the quasi-2D scattering amplitude /(Q^^) 
can be proved in the same approach. The definition of 

J(Q2D) ig 



/(2D) (c' ^ c) =-2^2 /^^^^(O) (^) I ^^^(^) |^( + ) (^^^ 



(D4) 



Using Eqs. ((M)) and ([29)), we immediately get Eq. ([39 



Appendix E: 2D scattering amplitude without SOC 

In this appendix we derive the 2D scattering amplitude 
of two distinguishable atoms in systems without SOC 
13811 . In this case we have the incident wave function 



^i°^ (P) - ^e^'-' 
Zn 



(El) 



with scattering energy e = fc^ . As shown in Appendix A 
and Ref. [33 , the scattered wave function ^^ ' (p) and 
the 2D scattering amplitude /'■^^■' (e) are given by 



^i^^ (P) = €^ (P) + Agnie; p, 0); (p > p.) (E2) 



and 



n-=0 



/(2D) (,) ^ ^^A. 



(E3) 
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Here, we have 



9r{£;p,p') ^ ±{p 



Ho-^ (P - P') 



^^^0 {-iV^p) , 



(E4) 



where iCg is the modified Bessel function of the second 
type. The coefficient A can be determined by considering 
that ^ (p) oc Inp— In din the region p* ^ p ^ 1/v^j with 
d the characteristic length. Using this condition and the 
expression ()E4|) of 5fl;(e;p, 0), we get the coefficient A 



and then the 2D scattering amphtude /q 



(2D). 



fr' (e) 



H'^^^ (0) 



2^ (hid 



c 



m^ 



(E5) 



where C = 0.5772... is the Euler gamma number. 

The expression (jESP clearly describes the character of 
the 2D low-energy scattering. It shows that the low- 
energy scattering amplitude /q (e) is totally deter- 
mined by the parameter d. In particular, /q (e) log- 
arithmically decays to zero [38| in the limit £ — >■ 0, and 
achieves the maximum when e = 4 exp (— 2C) /d^. 



Appendix F: Renormalization of the effective 2D 
interaction 



In this appendix we prove Eqs. (f53l [55 |) and (fSSj) for the 
renormalization of the effective 2D interaction potentials 
Vo and Vb- For the convenience of our calculation, in 
this appendix we do not work in the p-representation 
as before. We use the Dirac vector |)j_ to describe the 
quantum state of the 2D spatial relative motion of the 
two atoms, and \)t for the total quantum state for both 
the relative motion and spin of the two atoms. Namely, 
the total quantum \iP)t and the spinor wave function 
\ip (p)) defined in Eq. (1) is related as 



Wi 



dp\i^{p))\p) 



(Fl) 



where |p)_L is the eigen-state of the atomic relative coor- 
dinate. 

We first consider the single-channel potential Vo de- 
fined in Eq. (j52p . To obtain the correct expression for 
the parameter g, we should calculate the 2D scattering 



amplitude given by Vo- This scattering amplitude can be 
obtained from the Lippmman-Schwinger equation for the 
2-body T-operator To with respect to Vo'- 



to (E) ^Vo + Vo9 (E) Vo (E) , 
where the Green's operator g (77) is defined as 

1 



Hv) 



i] + iQ+-H, 



(2D)- 



(F2) 



(F3) 



and satisfies g {r], p, p') =j_ (p\g (77) \p')j_ with g (77, p, p') 
defined in Eq. P^ . Eq. (jF2p can be solved directly 
with the first-quantization form of Vo- We find that the 
on-shell element of To (E) has a separable expression 

T{i^[':^\fo{e,)\4'^)T^Vic')Uic), (F4) 

where \ipc )t is defined in Eq. ([6]) and we have Ec ~ 
£c' ■ The functions V (c') and U (c) can be obtained from 
Eq. (jF2p . The 2D scattering amplitude given by Vo is 



fo{c' ^ c) ^ -2^^ T{i^f\fo{e,)\^'^^> 



T- 



(F5) 



As shown in Sec. IV, the 2D scattering amplitude given 
by Vo should be the same as the quasi-2D scattering am- 
plitude, and then we should have 



/o (C' ^ C) = /(Q^I^) (C' ^ C) 



(F6) 



with /(Q2D) (jQgj^f;(j ijj Eq. (pg|) . This requirement di- 
rectly leads to Eq. ([55]) . 

Now we consider the two-channel model Vf, and prove 
the renormalization relations ([55]) and (|56p of Vb- To this 
end, we should calculate the bound state |^^*^ )t, which 
satisfies 



I* 



(eff)\ 



g{E,)Vb\^t''^)T. 



(F7) 



Equation (jF7p can also be solved directly via the first- 
quantization form of Vb- We can derive the energy Eb 
as well as the components of \^i^ )t in both the open 
and the close channels. As shown in Sec. IV, the en- 
ergy Eb of \W^^ )t should be the same as that of the 
quasi-2D bound state |^b)T, and then satisfies Eq. pB]) . 
The open-channel probability of |^|,'' )t is the same as 
the transverse-ground-state probability of \^b)Ti which 
is given in Eq. (|4ip . These requirements directly lead to 
Eqs. dnSl) and ^. 
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